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dom decoupling schemes, applicable to quantum dynamical control of arbitrary finite- 
dimensional composite systems. The basic idea is to design the control propagator according 

T 1 

^ ' to a random rather than deterministic path on a group. We characterize the performance of 

O ■ random decoupling protocols, and identify control scenarios where they can significantly 

weaken time scale requirements as compared to cyclic counterparts. Implications for reli- 
able quantum computation are discussed. 
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We present a general control-theoretic framework for constructing and analyzing ran- 
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5^. Dynamical decoupling offers a versatile control toolbox for quantum dynamical engineering 

in both traditional settings like high-resolution spectroscopy [ I ] and quantum information sci- 



I : e„ce [ ). DeccpUng .che^e. opera.e .bjec.ing .he ,.ge. .,.e. .o a series of open-loop 

control transformations, in such a way that the net evolution is coherently modified to a desired 
one [3]. This combines intrinsic design simplicity with the ability to avoid auxiliary memory and 
measurement resources, while additionally enabling straightforward integration with other pas- 
^ I sive [ ] or active [ ] quantum control techniques. Applications of decoupling range from the 
■ removal of undesired couplings in interacting quantum subsystems to active decoherence control 
and symmetrization in open quantum systems [6]. In particular, the use of decoupling meth- 
ods in conjunction with procedures for universal control [7] provides a route to noise-suppressed 
quantum computation based solely on unitary means. Remarkably, recent advances support the 
potential for highly fault-tolerant control architectures [ , ]. 

So far, general formulations of the decoupling problem have been restricted to deterministic 
control actions. In the simplest, so-called bang-bang setting, where the latter are instantaneous 
rotations drawn from a group Q, decoupling according to Q is enforced by cycling the control 
propagator through all group elements, translating into pulse sequences with minimal length Tc 
proportional to the size of ^ [ ]. This suffers from two main drawbacks. Because averaging re- 
quires traversing all of ^ in a suitable sense, decoupling becomes very inefficient for large groups, 
leading to unrealistically high control rates if the interactions to be removed have a short cor- 
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relation time r^. Furthermore, it is not clear how to handle interactions which are themselves 
fluctuating on time scales short compared to the averaging period T^. These limitations severely 
constrain the practicality of decoupling as a strategy for decoherence suppression in open systems. 

In this Letter, we propose to overcome the above limitations by introducing a framework for 
random dynamical decoupling. Physically, our approach takes inspiration from a naturally occur- 
ring instance of a random decoupling process; that is, the self- averaging of intermolecular inter- 
actions in gases and isotropic liquids due to random translational and re-orientational motions [i ]. 
This intuition is cast in control-theoretic language by requesting that the control propagator follows 
a random but known path on ^ [10]. We show how random decoupling may be used to achieve a 
desired coherent averaging and obtain a bound on worst-case performance. By comparing to or- 
dinary cyclic schemes, we find that in the presence of rapidly fluctuating interactions and/or large 
control groups, randomized design may prove superior. From the point of view of decoherence 
suppression, this not only establishes in general the counterintuitive possibility to actively cancel 
noise using randomness, but it also opens new prospects for significantly mitigating time-scale 
requirements in a wide class of control systems. 

Random decoupling setting. Let 5 be a quantum system with state space Ti^, AimiJ-Ls) = 
d < oo, evolving under an arbitrary, possibly time-dependent drift Hamiltonian Ho(t). Without 
loss of generality, we assume Ho{t) to be traceless for all t. We begin by constructing a random 
decoupling protocol for effectively switching off the evolution due to Ho(t), under the assumption 
of perfect, unbounded control. Let the available control generate a discrete or continuous compact 
group Q, acting on Hs via a faithful, unitary, projective representation /i, fi(g) = g for g E G, 
/i(^) = ^. A random decoupler uses control in Q in two ways: first, to establish a logical frame 
that is related to the physical one (where i?o(^) is specified) by an element of Q; second, to rotate 
the system according to Q randomly over time, by following a random control path Uc{t). Thus, 
both the past control operations and the times at which they are applied are known, but the future 
control path is random. 

The essence of the random decoupling approach is to directly depict the evolution of the system 
in the logical frame that continuously follows the applied control. Let psit) = U{t)ps{Q)W{t) 
describe the state of the system in the physical frame, evolving under the action of both the internal 
Hamiltonian and the controller, and let psit) = Ul{t)ps{t)Uc{t) denote the corresponding logical 
state, with ps (0) = P5 (0) . Then the evolution in the logical frame is fully specified by a propagator 

{h=i) 

tj{t) = Ul{t)U{t) = Texp duH{u) 

where H{t) = U}.{t)Ho{t)Uc{t). Under the usual cyclicity assumption of deterministic decou- 
pling, Uc{t + Tc) = Uc{t) for Tc > 0, the physical and logical frames stroboscopically coincide at 
times tj^ = NTc, G N. By contrast, random decoupling is intrinsically acyclic, and the control 
path almost never returns the system to the physical frame. However, the available information 
about the past control trajectory may be exploited to bring the state of the system back to the 
physical frame if desired. 

Error bounds. To determine whether and how well random decoupling succeeds at suppressing 
the dynamics due to Ho{t) it is necessary to compare the evolution under the propagator U (T) over 
a time interval T to the identity evolution, up to a global phase. A natural measure is provided 
by the error probability for an arbitrary pure initial state Ps = |^)(^| of S. With respect to the 
random nature of the control path, the a-pnori error probability can be expressed as an expectation 

eriPs) = E {tr^ {P^ps{T)) } = E [tvs (^P^U{T)PsU{Ty) } , (1) 
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where Pg = t s — \ip) {ip\ is the orthogonal complement of P5 and E denotes ensemble average. 
Then a worst-case pure state error probability may be defined as 

eT = Maxp,,{er(P5)} • (2) 

A quantitative bound for er is contained in the following: 

Theorem 1. Suppose that (i) Q acts irreducibly on Hs- (n) Udf) is uniformly random for 
each t. {iii) For any t,s > 0, Ucif) and Uc(t + s) are independent for s > At. {iv) ||ifo(^)||2 is 
uniformly bounded in time by k > 0. Then 

eT = {TAt k'^) for TAt A;^ < 1 . (3) 

Here, \\A\\2 = Max |eig(\/y4M)|, and uniformly random is intended relative to the invariant 
Haar measure on Q, normalized such that i'g{Q) = 1 ]. While a rigorous proof of the 
above Theorem is rather lengthy [ ], an outline of the underlying strategy suffices for gaining 
physical insight. The key step is to realize that, in each of the integrals involved in the Dyson series 
expansion of the time-ordered exponentials defining U (T) and f/^(T) in Eq. (1), the independence 
assumption (Hi) effectively partitions the integration domain in two separate regions: a volume 
14^1 (At), where none of the integration variables is more than At away from all the remaining 
ones; and the complement W2{At), where this condition is violated by at least one variable. The 
expectation relative to such a variable may be taken separately, leading, under the uniformity 
assumption (ii), to a contribution of the form 

E{Ul{t)Ho{t)U,{t)} = [diygg^Ho{t)g, t G (0, At) . 

Jg 

Notice, as the result of such an ensemble average, the appearance of the same dynamical Q- 
symmetrization which, in standard deterministic schemes, is achieved through the time average 
over a cycle [3, 13]. In particular, the irreducibility assumption (i) implies maximal projection in 
the set of scalars. That is, for X traceless, 

j^dugg^Xg=^^ts = ^- (4) 

As a consequence, all terms originating from 14^2 (At) vanish, and the desired upper bound to 
^t{Ps) may be determined by estimating the volume of lyi(At). The irreducibility assumption 
can of course be weakened. As it turns out, the final result (3) for has a simple intuitive 
explanation, which we defer until after we describe the corresponding error bound for deterministic 
schemes. 

From an implementation perspective, one may distinguish two main scenarios, depending on 
whether the decoupler is specified by a continuous or discrete control group Q. In the former 
case, the decoupling time scale At is defined by the independence requirement between Uc{t) and 
Uc{t + s), condition [iii] entering as a design constraint. Note that bounded-strength controls 
might suffice as long as At is finite. If Q is discrete, the required random walk of Uc{t) may 
be enforced through a sequence of equally spaced bang-bang pulses randomly drawn from Q. 
In this case, the independence requirement is automatically satisfied by identifying At with the 
separation between consecutive kicks. Either way, it is important to stress that random decoupling 
(unlike deterministic decoupling) places no restriction on the temporal behavior of Ho(t), only on 
its maximum eigenvalue. 
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Random decoherence suppression. The above formalism can be extended to the suppression of 
noise effects arising from the coupling between the target system S and an uncontrollable quantum 
environment E. Let the total drift Hamiltonian be expressed in the form HQ{t) = I5 ® He + 
J2a Ja{t) ® Ba, whcrc He accounts for the (typically unknown) evolution of E and the internal 
evolution of S is included among the interaction operators, with tr( Ja(t)) = for all t. The action 
of the decoupler is understood as Uc{t) ® 1e. Physically, it is meaningful to define a pure-state 
error probability that depends only on the reduced state of 5 in the logical frame. That is, ps{T) 
in Eq. (1) is now calculated as ps{T) = tTE{U{T)psE{0)W (T)}, Pse{0) = Pse{0) being the 
joint initial state and, as before, the logical propagator U{t) describing the combined evolution in 
a frame that explicitly removes the control field. By purifying the environment, we can assume 
that Pse{Q) = -Ps ® -Pb, both P5 and Pe being one-dimensional projectors. The derivation of a 
bound for eriPs) be formally carried out following the same steps as in the uncoupled case. 
It suffices to observe that Eq. (1) is equivalent to 

eriPs) = E [trs,E [P^ ® tEU'{T)Ps ® PeU'ITY) } , 

with the propagator 

U'{t) = Ul{t)Ul{t)U{t) = Texp duH'iu)"^ 

describing the evolution in a frame where both the applied control and the environment dynamics 
UEit) = exp{-iHEt) are explicitly removed, and H'{t) = Y,a Ul{t)Ja{t)Uc{t) ® Ba. We thus 
have the following: 

Theorem 2. Let Q act irreducibly on Tis and satisfy the same uniformity and independence 
assumptions as in Theorem l.lf\ \ Yla Jaif) ® -Ba| I2 is uniformly bounded in time by \> then 

eT = {TAt A^) for TAt < 1 . (5) 

Formally, A is a measure of the overall noise strength as defined in the context of quantum error 
correction theory [14]. As pointed out in this reference, caution is required in treating infinite- 
dimensional environments. Physically, 1/A = Tc is of the order of the shortest correlation time 
scale present in the interaction to be removed. While the latter provides the relevant time scale to 
consider in the absence of additional information about the environment's initial state, power spec- 
trum, and internal dynamics, such properties may critically impact the decoupling performance in 
actual applications [15]. Thus, lower error bounds tend to be fairly example specific. 

According to the above Theorems, er can in principle be made arbitrarily small by appropriate 
control design, implying the possibility to arbitrarily suppress on average the unwanted evolution 
in the logical frame. This is especially surprising for decoherence suppression considering that, 
in the physical frame, the applied random field appears to be in general a source of decoherence. 
It is worth noting that the possibility to exploit randomization was considered earlier for specific 
decoupling problems. Preservation of coherence of a lossy radiation mode via the random mod- 
ulation of a system parameter was established in [16]. More recently, a randomized refocusing 
algorithm was proposed in [17] in the context of efficient simulation of quantum computation 
starting from few-body Hamiltonians on n qubits. While revisiting such specific situations in the 
light of the present analysis is interesting in itself, our main goal in what follows is to continue 
developing a model-independent formulation of random decoupling in general control-theoretic 
terms. 
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Comparison with cyclic decoupling. In order to assess the performance and usefulness of ran- 
dom decoupling schemes, a comparative error bound for deterministic decoupling is needed. We 
focus on the standard situation where the drift Hamiltonian Hq is time-independent, and decou- 
pling is accomplished by cyclic averaging over a finite group of order \Q\ > 1. Apart from the 
redundant ensemble expectation, Eqs. (l)-(2) still define a valid worst-case pure state error prob- 
ability. The deterministic counterpart to Theorem 1 is then the following: 

Theorem 3. Suppose that (i) Q acts irreducibly on Hs- {ii) Ucif) is assigned according to 
a cyclic path over Q, with Udf) = gj for t G [j At, j + lAt), j = 0, . . . , |^| — 1, At > 0, and 
Tc = 1^1 At. (m) IliJolb ^'■5' bounded by k > 0, with kTc < 1. Then 

er = O {{TT.k^f'^ for TT,k^ <^ I - kT, . (6) 

The proof follows from a direct evaluation of the logical propagator U (T) using average Hamil- 
tonian theory [ ], 

where H is computed from the Magnus expansion under the averaging and convergence condi- 
tions, H^^^ = and kTc < 1, respectively [12]. We now provide an intuitive justification to the 
error bounds we found. 

Write R = TAtk'^ = (fcAt)^(T/At). For the random method, each control step can accu- 
mulate an error amplitude of up to kAt. Randomizing the decoupler has the net effect that the 
amplitudes add up probabilistically. Therefore, over an evolution time T, the total error probabil- 
ity is bounded by the number T/ At of such intervals, times the error probability [kAtY of each 
step. Notice that the bound of Theorem 1 is indeed = 0{R). 

For the cyclic method using \Q\ steps of duration At in each cycle, the dominant errors are due 

to H^^\ That is, they arise from non-commuting contributions associated with pairs of intervals 
in a cycle. Thus, for each cycle the error amplitude is bounded by |^p(/cAt)^, and a total time 
T contains T/(|^|At) such cycles. If, as assumed, each cycle is identical and the interaction is 
constant, the total error amplitude is bounded by the sum, yielding \Q\R. By squaring and using 
that 1^1 At = Tc, the bound of Theorem 3 emerges, = Od^pi?^). 

The above analysis shows that the worst-case errors of the two procedures compare as follows: 

= 0{R) vs = O {{\g\^R)R) , 

the quantity \Q\'^R becoming a relevant figure of merit for performance. Thus, cyclic decoupling 
tends to perform better if any time dependence or fluctuations in the interactions to be removed 
have time scale longer than = |^|At and, in addition, \g\'^R ^ 1. Superior performance of 
random decoupling is expected instead in situations where the effective correlation or fluctuations 
have time scales large compared to At but short compared to |^| At; or, alternatively, ^ 1. 

Generalizations and applications. The above results lend themselves to a number of general- 
izations. The extension to reducible group actions (hence selective decoupling) is conceptually 
straightforward. Procedures for universal decoupled control may be designed similarly to [7], by 
randomly modulating the applied control Hamiltonians to compensate for the decoupler action 
if necessary. This paves the way to schemes for randomly controlled noise-suppressed universal 
quantum computation. In addition, one may envisage a variety of hybrid control schemes where 
deterministic and random operations are simultaneously exploited. At least two options are worth 
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considering. First, one may randomize the decouplers. If multiple decouplers are available to 
effect a desired averaging, which one to apply may be picked at random at every cycle. Or, with 
a single decoupler, one may randomize the cycles, by randomly choosing which path to follow 
to traverse Q. While a clever concatenation of deterministic and random protocols could merge 
advantageous features from both methods, quantitative error estimates as well as studies of the 
typical performance in specific situations will be reported elsewhere. 

We anticipate that randomization might offer substantial benefits whenever a large number of 
control time-slots is involved. An extreme example is maximal decoupling in n arbitrarily cou- 
pled qubits, d = 2". Deterministic group-based schemes require averaging over the Pauli error 
basis {]1, 0"^., (Ty, a^}*^", with \Q\ = d? = [3]. For fixed control parameters T, At such that 
R <^ \, the condition \G\'^R <^ 1 becomes exponentially harder to meet as n increases. Equiva- 
lently, for a fixed tolerable error e^, an interval At that shrinks exponentially with n is needed to 
compensate \Q\'^ in this case. A randomized implementation of Pauli decoupling is indeed at the 
heart of the simulation algorithm mentioned above [17]. In addition, the recently proposed Pauli- 
Random-Error-Correction method for coherent errors [IS] may also be understood as an ingenious 
application of the present control framework, random Pauli rotations being repeatedly applied to 
average static imperfections, and permutations of the original logic gates ensuring the intended 
decoupled control. While cyclic schemes with quadratic complexity [19] are known for bilinearly 
coupled qubits as assumed in [18], randomized schemes may still be attractive for large n and/or 
time-varying couplings. In the same spirit, the cancellation of rapidly fluctuating dynamical im- 
perfections reported in [20] may be suggestively reinterpreted as a random self-decoupling effect. 
Lastly, efficiency improvements are to be expected from decoupling according to the symmetric 
group Sn acting on n qubits, which otherwise involves factorial overheads, and is relevant to the 
synthesis of collective noise [^ ' ]• 

Conclusion. We introduced an approach to dynamical decoupling that relies on random control 
design. Beside being interesting per se as a largely unexplored setting for coherent and error 
control, random dynamical decoupling carries the potential for faster convergence and relaxed 
timing constraints compared to deterministic counterparts in relevant situations. While additional 
work is needed to expand the present analysis, we believe that our results add to the significance 
of decoupling methods as a control-theoretic tool and allow a step forward toward making them a 
practical error control strategy in quantum information science. 

L. V. acknowledges support from LANL during the early stages of this work. E. K. was sup- 
ported by the US NSA. We thank Howard Bamum for discussions. Contributions to this work by 
NIST, an agency of the US government, are not subject to copyright laws. 
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Appendix 



In this Appendix, we supply the proofs of Theorem 1 and 3. Theorem 2's proof can be obtained 
from that of Theorem 1 by calculating the relevant worst-case error probability in the appropriate 
doubly-rotating frame. 

Proof of Theorem 1. Let Ps be an arbitrary pure state of S. The logical propagator U (T) may 
be expressed as follows: 

oo 

t/(T) = ^/„(T), (7) 

n=0 

where 



/„(T) = i-tr / dur,... du.HM . . . Hiu,) , (8) 

and similarly for U (r)^ Thus, we need to calculate 

eTiPs) = ^Urs( f; P^IniT)PsImiTA \ = f; E{trs{P^IniT)PsUTy)} . (9) 

^. \n,m=0 / J n,m=0 

The contributions with n = or m = vanish because of Pg and Ps cancel each other upon 
exploiting the cyclicity of the trace. By noticing that exiPs) > hence €t{Ps) = l^riPs)]^ 

eT{Ps)< Yl \^{i^s{PsUT)PslUTy)}\ . (10) 

n,m>l 

Under the assumption of sufficiently smooth behavior, the expectation may be moved under the 
integral. Fix a pair of integers n,m> 1, then the relevant contribution is 

/ du^...dUndh...dt^¥.{p^H{un)...H{u^)PsH{h)...H{t^)] . 

JO<...<Un<T;0<...<tm<T ^ ^ 

(11) 
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Let W^'^''^\At) denote the set of points (ui, . . . , u„, ti, . . . , tm) satisfying that and t£ are each 
time-ordered and no ug or is further away than At from the rest, and let l^g"'"^^ (At) denote the 
remaining integration volume in Eq. (11). Because, within W2'^''^\At), at least one of the inte- 
grating variables is more than At away from all the other variables, the independence assumption 
(Hi) allows the expectation relative to such a variable to be taken separately. By the uniformity 
assumption on Uc{t) for all t, and by the tracelessness assumption on Ho{t) for all t, such an expec- 
tation vanishes. Therefore, W^'^'"''\At) is the only subset of points contributing to the expectation 
in Eq. (11). Let rfw^"'™) denote the corresponding integration measure. Then 

eAPs) < E / ^^^"'"^ h {^'' {^sHM ...Ps... H{t,n)) } 

n,m>l ■^M^i"'™ 



(12) 



where in the second step Jensen's inequality has been used. By noticing that the argument of the 
trace is a rank-1 operator, one can simplify 

eriPs) < E / dw^'''"'^E\\\P^H{un)...Ps...H{tm)\\2} 

< y / dM;("'"^)E{||P^||2 llPsllsA;"^™} 

< J2 voi(iyi^"''"^)r+'" , (13) 

n,m>l 

where the uniform bound k for ifo(^) has been used, and Yo\{Wi^'"^^) is the volume of Wi"''"^\ 
Note that the dependence upon the initial state Ps has disappeared at this point. 

The above volume may be estimated through the following combinatorial argument. First, 
notice that given the two ordered lists < ui < . . . < Un < T , < ti < . . . < tm < T , there 
are ("^'") different merged orderings. Fix a particular one. Then each element needs to be either 
within At of the next one or of the previous one. Make a choice for the odd-numbered elements, 
the first element being labeled 1. There are at most 2^*^"+™^/^^ such choices. For each of them 
the contribution to the volume may be bounded by ordering the even-numbered elements, then by 
inserting the odd ones, ignoring the ordering constraint now. Finally, 

voi(iy^™^) < (^ + '^)2r(«+™)/2i^ 



m J {[{n + m)/2\)\ 

^ 2l{n+m)/2']rpl{n+m)/2] i^2/\t)^^'^~^'^^^'^\ (14) 

where the inequalities < 2"+™"^ (for n + m > 2), and [{n + m) /2J ! > 2L("+'")/2J-i have 

been exploited. 

The last step is to sum over n, m: 

oo 

Maxp^ieriPs)} = er < ^ 2r("+"^)/2lTL("+'")/2J (2At) . (15) 

n,m=l 
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This may be done by considering separately the four partial sums where both n and m have the 
same (even or odd) parity, or they have opposite (even-odd or odd-even) parity, respectively, and 
by evaluating the [J, [] in Eq. (15) accordingly. Straightforward calculations yield 

er < (4TAtfc^) (1,4^^,^2)2 = OiTAtk') (16) 

for values of TAtk'^ <^ 1, as quoted in Eq. (3). □ 

Proof of Theorem 3. The logical propagator U (T) may be expressed in terms of the average 
Hamiltonian H as 

00 _ 

U{T) = e-'^^ = V -i-tHTT , (17) 
^-^ n\ 

n=0 

and similarly for U (T)^ Then the desired pure-state error probability is bounded by 



00 



eriPs) < E {P^i-tHTrPs{+tHTr)\ , (18) 



n.m=l 



where the orthogonality of the projectors Ps, Pg has been used to remove the terms with n = 
and/or m = 0. By observing that the argument of the trace is a rank-1 operator this gives 

eT{Ps)< V ^^||P5^||2||P||2||:ffT||^+'"< (l-ell^ll^^) . (19) 

n,m=l 

Assuming that first-order averaging and convergence conditions for the Magnus series are fulfilled, 
one has H^'^^ = and kT^ < 1, implying 



00 



1 - kT^ 
e=i ^ 



In the limit of sufficiently short time, k'^TT^ < 1 — kT^, | I2T < 1 and by using the inequality 

11 — e^^l < |x|/|l — x| for X < 1 one obtains 



Maxp,{eHP,)} = er< = O ((| ^T)^) = O ((TT.A:^)^) (20) 

for TT^P < 1 - kTc, as quoted in Eq. (6). □ 



